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Abstract — We develop a new statistical test for the wavelet 
power spectrum. We design it with purpose of testing signals 
which intrinsic variability displays in a Fourier domain a red- 
noise component described by a single, broken or doubly- 
broken power-law model. We formulate our methodology as 
straightforwardly applicable to astronomical X-ray light curves 
and aimed at judging the significance level for detected quasi- 
periodic oscillations (QPOs). Our test is based on a comparison 
of wavelet coefficients derived for the source signal with these 
obtained from the averaged wavelet decomposition of simulated 
signal which preserves the same broad-band model of variability 
as displayed by X-ray source. We perform a test for statistically 
significant QPO detection in XTE J1550-564 microquasar and 
active galaxy of RE J1034+396 confirming these results in the 
wavelet domain with our method. In addition, we argue on the 
usefulness of our new algorithm for general class of signals 
displaying l//"-type variability. 

Index Terms — wavelet transforms, signal analysis 



I. Introduction 

Fourier transform has become a very popular tool of time- 
series analysis within last fifty years in almost all areas of 
research. Applied to test signal frequency content, it was 
quickly recognized as a helpful method for investigation 
of periodic variability |T). It gained a particular interest 
in a wide field of astronomical research regarding different 
classes of objects with variable emission. In particular, time- 
series provided by recent high-energy detectors onboard X- 
ray orbiting NASA/ESA satellites 0-113] have been found 
to trace remarkably rapid variability coming from neutron 
star and black-hole systems |4|-|6|. Since the majority of X- 
ray emission takes its origin in accretion process of matter 
and gas onto central object Q-JSl, therefore the observed 
signal variability is a direct manifestation of physical processes 
taking place in a very strong gravitational field JO]. 

The accretion onto black-holes is one of the most energetic 
process in the Universe LlOj . Energy spectra for the corre- 
sponding objects uncover both thermal and non-thermal nature 
of emission (e.g. 1 1 1 1). The analysis of Fourier power spectra 
(also referred to as Power Spectral Density; PSD) of X-ray 
light curves x{t) has revealed that the overall shapes are very 
often well described, in general, by a power-law model, i.e. 
P{f) oc 1//" where a G denotes a slope ([31, fT2l-fT6l). In 
case of a > one claims about signals dominated by intrinsic 
variability of red-noise type. The red-noise character of source 
variability is generally rarely met in nature |17|. Nevertheless, 
it seems that black-hole systems emitting in X-rays enlarge 
the sample significantly. 
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To date the majority of efforts devoted to our curiosity of 
understanding the X-ray variability of astronomical objects 
were based on the systematic study of PSD shapes and looking 
for occasionally observed quasi-periodic oscillations (QPOs; 
f5\, ifTsl -lfTo!). Since Fourier PSD contains no information 
about time evolution of detected periodicities, the need of 
use of the time-frequency techniques emerged. The application 
of Short-Term Fourier Transform (STFT; e.g. f20l) provided 
a good solution in this domain and has been found helpful 
in localizing QPOs (e.g. ll2Tl - ll22l '). However, regardless the 
frequency of dectected QPOs, STFT keeps the same time- 
frequency resolution as dictated by a length of a sliding 
window. In most cases it prevents a detection of oscillations 
lasting shorter than the window span. On the contrary, a 
wavelet analysis occurred to be more successful where the 
time-frequency resolution was scale (frequency) dependant. 
Within recent fifteen years it has attracted attention of many 
researchers Il23l - ll24l . 

The wavelet power spectrum (scalogram) for discrete evenly 
sampled X-ray signal x„ (n = 1,...,N; T = NAt) given as 
count-rate [cts s~^] can be defined as the normalized square 
of the modulus of the wavelet transform: 



W = f |w„(a,„) 



(1) 



where ^ denotes a normalization factor and w a discrete form 
of the continuous wavelet transform in a function of two 
parameters: scale a and localized time index n ||24| - ||26I : 
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where the discrete Fourier transform of signal Xn is given by: 
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= '^{xn~ x)e 
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and j denotes a frequency index, 

^ ^ f j/iNAt) for J < N/2 
\~j/{NAt) forj>7V/2 



(3) 



(4) 



and x a mean value of Xn ll26l . Following ll27l one can 
assume normalization factor of w„(am) to be ^ = 2Atx~^ 
which provides that integrated wavelet spectrum over time, 
also referred to as a global wavelet power spectrum 
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(5) 



will hold the units of (rms/mean)^ Hz ^, i.e. it gives the 
variance, relative to the mean, within a given frequency range 
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of integration ||27ll - l|28l . For the computation of ([TJ the use of 
a dyadic grid of scales, i.e. 

a™ = ao2'"^'", m = 0,...,M, (6) 

can be implemented |!26|, where the smallest scale aq = 2At 
corresponds to the reversed Nyquist frequency, and 

M = Am-^ \og^{NAta^^). (7) 

A grid resolution is given by Am. Since the search for 
periodicity and quasi-periodicity in X-ray light curves seems 
to be the main target, it is reasonable to assume the analyzing 
function to be Morlet wavelet, 

^(i)=^-l/V(2./ot)g-tV2^ (8) 

which oscillates due to a term oc e'*, where 27r/o parameter 
is set to 6. 

In order to popularize wavelet analysis as a new tool in 
research over time-series, an additional requirement in the 
form of significance test was strongly desired. Torrence & 
Compo |26| (TC98) were first who addressed that issue 
extensively providing the public with an excellent guide and 
software ready-to-use. In general, a map of wavelet power 
spectrum shows a distribution of spectral components in the 
time-scale (or time-frequency) plane. Single oscillation, of 
sufficiently high amplitude to be detected, marks itself in the 
map as a wavelet peak. To every peak a certain statistical 
significance can be assessed. We say that a wavelet peak is 
significant at assumed per cent of confidence when it is above 
a certain background spectrum. The latter can be defined by 
the mean Fourier power spectrum of analyzed time-series. 

TC98 in their work adopted the background power spec- 
trum, Pj, after [291 (see also Eq.(16) in |26|) which is valid 
for application as long as the analyzed signal is considered 
to be a realization of univariate lag-1 autoregressive (ARl) 
process li30J -l.32l. xi = axi-i + zi, where a is lag-1 
autocorrelation, zi denotes a random variable drawn from the 
Gaussian distribution (with zero mean and variance ct^) and 
xo = 0. However, the PSD shape provided by 1291 is not a 
good representation of broad-band PSDs in case of considered 
here X-ray sources. This is because of at least two reasons: 
(i) not every X-ray light curve can be described by ARl 
process |[33l (although 1341 ), (ii) PSDs are dominated by the 
red-noise spectral components which shapes are mainly fitted 
with single, broken, doubly broken or even bending power-law 
models. 

TC98 showed that for a time-series modeled as the ARl 
process its local wavelet spectrum follows the mean Fourier 
spectrum given by their Eq.(16) ||26l . By the "local wavelet 
spectrum" they define a vertical slice (along the frequency 
axis) through a wavelet map. Empirical justification is done 
on the way of Monte Carlo simulations which were performed 
considering only one local wavelet power spectrum per sim- 
ulation and, in addition, selected for a particular moment of 
time (k = 256 of 512 points). The corresponding distribution 
of the local wavelet spectrum (at each localized time index n 
and scale &„), ought to have X2 distribution, 

^ Ip^xI (9) 



where X2 denotes the distribution with two degrees of 
freedom and an arrow means "is distributed as". The 

above implication is performed based on the assumption that 
a random variable xi is normally distributed. 

A noteworthy work in this domain has been recently done 
by Zhongfu Ge [35 1 who generalized TC98's significance test. 
Ge left the assumption of signal to be the realization of ARl 
process and proposed Gaussian White Noise process to serve 
as the null hypothesis. 

What TC98 and Ge did not show is that the local wavelet 
power spectra, Pj, at nearby times n are correlated due to 
variable length of the wavelet function, and the correlation is 
stronger at low frequencies as shown by l36ll . This introduces 
departures off distribution in any finite data set. Thus a 
more correct procedure is requked to take this effect into 
account. 

Maraun, Kurths & Holschneider fyj] proposed another test 
for wavelet significance, namely, areawise test which aims 
at comparing the wavelet peak size to the expected time- 
frequency uncertainty. They noticed that for any process with 
unknown distribution a number of wavelet peaks forming in 
the wavelet map different shapes (circles, ovals, patches, etc.) 
become difficult in interpretation. The extracted picture is 
messy. Some peaks are due to red-noise process of the source 
intrinsic variability and the others are due to white-noise. As 
shown by them, the areawise test may reduce a number of 
significant peaks even by 90 per cent. 

To judge the reality of each wavelet peak it seems best 
to set the significance level highly enough in order to reject 
spurious results. However, not always this is a right way, for 
instance, in case of X-ray sources of rather low count-rate 
(low S/N ratio) and strong mixing between red- and white- 
noise components in high-frequency region. On the other hand, 
this also does not mean that the analysis does not contain 
any scientific results. For example, the same situation took 
place for COBE detections of microwave background peaks 
||381 . Most of them were spurious but the fluctuations were 
statistically different from expectations of pure noise f39l. 

None of this method can be applied as a fully trustful 
test working well for astronomical X-ray signals as long as 
their assumptions rely on the analyses that use impropriete 
background spectra. In this paper we address that issue and 
propose a new significance test. 

II. Wavelet significance for red-noise signals 

OF POWER-LAW FORM 

We design a method for determination of significance of 
wavelet spectral features (in short, wavelet peaks) calculated 
for any X-ray light curve which intrinsic variability is a real- 
ization of a sum of red- and white-noise (Poisson) processes. 
The main idea standing for that is to compare the distribution 
of wavelet power between real data and simulated signal. Here, 
by the simulated light curve (also referred to as TK light curve) 
we will understand a time-series which displays the same 
profile of red-noise variability as that one found for X-ray 
data. Before providing the reader with ready-to-use algorithm 
(Sect. III-CI ), some essential information on the background 
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assumptions need to be clarified. We provide them within the 
following two subsections (Sect. UKAI and HTBI i. 

A. Simulations of 1//" process 

Let's assume that x{t) describes a part of stochastic con- 
tinuous process X{t) generated by the physical system of 
unknown nature which produces the variability. Based solely 
on x{t) a reproduction of the overall information and prop- 
erties of the system is very difficult. Recorded data can 
facilitate a determination of a model which would reproduce 
the system best. The observed light curve of X-ray emission is 
the realization of the physical process(es) for which calculated 
Fourier power spectrum (PSD) reveals the variability of power- 
law l//"-form (see references in Section Therefore, we 
are able to approximate one of the system's parameters (i.e. 
its variability) by a power-law model which, as we believe, 
constitutes the best PSD description of the underlying process. 

Aperiodic broad-band variability as displayed in X-ray light 
curves can be thought, in the first approximation, as the real- 
ization of a linear stochastic process, weakly non-stationary. 
Timmer & Konig |40| proposed an algorithm (hereafter also 
referred to as TK) for generating linear aperiodic signals, l{t), 
that exhibit P{f) oc 1//" power-law spectrum, 

lit) cx J2 /P(7) cos pvr/ - , (10) 
/ 

where the phase is randomized (/){/) G [0,27r]. Since the 
amplitudes are taken as y/Pjj), [40] noticed that in order 
to create power-law signal it is essential to randomize also the 
amplitude. 

A light curve obtained in this manner displays 1//"- 
shape, however as argued by li33l . it does not preserve 
fundamental statistical properties of X-ray time-series for 
accreting black-hole systems. The latter have been found to 
be rather non-linear showing rms-flux relation at all time- 
scales 1331 . 1411 (rms is defined as a square root of variance 
calculated for a signal). According to [33 1, formally non-linear 
type of variability for any linear input light curve can be 
obtained by taking simple exponential transformation of ( fTOl i. 
i.e. Xexp{t) = exp[/(t)]. To make the picture complete every 
simulated light curve should include white noise. Following 
II25I it can be expressed as: 

x{t) ^ POlSDEY[Xcxp{t)At]/At (11) 

where x{t) denotes simulated signal with Poisson noise. At 
is a light curve resolution (bin time) and POISDEV represents 
a subroutine generating a random number with Poisson distri- 
bution Il42l . 

B. Selection of best 1//" model for PSD of X-ray source 

A description of Fourier PSD of X-ray object under inves- 
tigation can be obtained in the process of numerical fitting 
of power-law model. As reviewed by 15], a small but sub- 
stantial subset of X-ray black-hole (and neutron star) systems 
apart from exhibiting l//"-like broad-band variability, also 
show quasi-periodic oscillations (QPOs). Hardly ever met 
as tiny spikes, they are much often resolved as broadened 



features characterized by the quality factor of Q = f/Af 
where A/ denotes the peak's full width at half maximum 
(FWHM) and / QPO central frequency. The broader peak 
the higher A/ thus the lower Q. So far, the majority of 
QPO peaks have been found to be surprisingly well described 
by Lorentzians, e.g. fT2l . lfT6l . Lorentzians are defined as 
L{f) oc Af/[f - fof + (A//2)2] with a centroid frequency 
/o and, if detected, they constitute an intergal component of 
the overall PSD model. In practice, fitting of power-law model 
as a background model may or may not include fitting of L{f) 
component simultaneously. In the latter case, PSD data points 
describing Lorentzian can be simply omitted and these two 
aforementioned approaches rely on the individual preferences 
and fitting experience. 

In this work our wish will be to generate simulated X-ray 
light curve which will own the same power-law background 
model of variability as fitted to the object's PSD. This step 
requires that the overall rms and the average count rate for 
simulated time-series must match the corresponding same 
values for X-ray signal. For Fourier PSD normalized to unit 
of (rms/mean)^ Hz^^ |27|, the root-mean-square (rms) can be 
obtained by integration of underlying PSD in a given band of 
frequencies l28l as follows: 

rms2 = = P{f )df « V P,Af (12) 

where 1/NAt < /i < /2 < 1/2 At. We argue that rms values, 
both obtained for simulated signal and X-ray light curve ought 
to be calculated in the same (/i, /2) range of frequencies what 
can guarantee the best match. The integration (fT2] | should 
omit PSD data points corresponding to QPO component(s), 
if present. 

C. A new algorithm for wavelet significance 

We define a new algorithm leading to determination of 
significance matrix for derived wavelet power spectrum by the 
following steps: 

> Two-dimensional matrix of wavelet power spectrum ([T]l 
for a time-series x„ is a set of (A/ + 1) x wavelet 
coefficients where M corresponds to (|7]i and N denotes 
a number of data points. Let us denote the results of 
computation of ([T]i for X-ray source as: 

Please notice that in practice instead of time-scale plane 
one often refers to wavelet spectrum displayed in time- 
frequency plane. If required, for Morlet wavelet, a trans- 
formation from scales to Fourier frequencies can be 
applied according to fj = (LOSom)"^ relation ll26l . 

> For a given discrete signal Xn calculate its Fourier power 
spectrum and find the best model which describes it 
(Section III-Bb . If the model for considered data set is 
known in advance (e.g. from the existing literature), use 
the best fitted parameters as the reference ones. 

> Generate one long simulated light curve (with the same 
At), best at least two or three orders of magnitude longer 
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than your signal under analysis. Use Timmer & Konig's 
method which allows to obtain simulated signals from 
Fourier power spectra of assumed shape (Section III-AI ). 
Take care to keep the rms variability at the same level as 
in the data (Section HTbI ). 

From the simulated TK light curve, select randomly K 
signals of the same time duration as in original data set, 
i.e. iVA<, and compute for each of them the wavelet 
power spectrum according to ([T]l. Shall us denote the 
resulting matrices as: 



and /: 



WJ 



where 



(14) 



(M+l)y.N 

In calculation of ( fT4l l for normalization factor ^ (see 
Sect, m assume every time a mean count rate x of the 
local TK signal i. 

Find the maximum value between all W"^ matrices: 



V 3 X where x ~ max(Vl^' 



(15) 



For every scale m of every calculate the histogram of 
wavelet power distribution from to x with a histogram 
resolution defined as /S.h = x/{M + 1). In consequence, 
the following matrices shall be created: 



{M+l)x(x/\h) 



(16) 



consisting of M + 1 histograms. 

Calculate a matrix containing averaged histograms for 
every scale m: 



(M+l)x(x/\h) 



1 ^ 

K ^ 



(17) 



This step is introduced in order to smooth out different 
distribution of wavelet power appearing in every for 
each scale. Please note the higher K the more smoothed 
averaged wavelet power histograms one obtains. 
For every of M + 1 histograms stored in H, determine 
the p-th quantile Xp. Namely, let H{m) be a continuous 
random variable, then for < p < 1, the p-th quantile is 
the Xp such that PT{H{m) < Xp) — p f44l. For instance, 
in deriving of the significance of wavelet peaks at the 95 
per cent confidence level one aims in determination of 
0.95-th quantile Xp in a function of scale m (frequency). 
In practice, a numerical integration of normalizeqj H{m) 
must be performed until an integral value exceeds 0.95. 
This determines quantile Xp with an error of A/i for each 
histogram: 



Xr, 



(18) 



(M+l)xl 



• Expand matrix Xp to have dimensions [M + 1) x by 
calculating a linear algebraic product of the matrices Xp 

'a term normalized means that an integral of H{m) must equal 1. 



Tj/bkg 



1 

E 

k=l 



Xp{i,k)I{k,j), (19) 



where / is a matrix of ones of 1 x iV dimension. Denote 
the result as a background wavelet matrix Vt^'^'^s 
• Finally, a matrix of (1 — p)100% significance, S, will 
identify significant wavelet peaks appearing in computed 
wavelet power spectrum W^'^'^ when the following con- 
dition will be met: 

^src 

'S'(M+l)x7V„b, = -^^Tbki > 1- (20) 

In proposed procedure a distribution of wavelet power at 
a given wavelet scale (frequency) is directly compared to the 
distribution of power obtained for simulated light curves. Since 
the Fourier PSD model and rms variability in simulated time- 
series are assumed to match the same quantities as for the 
X-ray source data, they serve as a good reference level to 
test any excess of wavelet power above it. In addition, by 
testing a wavelet power scale by scale, one takes into account 
an influence of correlation of wavelet power over a range of 
frequency and time due to varying length of a wavelet function 
(Sect. ID [36]). 

III. Application to X-ray observations 

In this section we provide the reader with two examples of 
application of our new algorithm to real observations of black- 
hole systems emitting in X-rays: microquasar XTE J1550- 
564 ESI and active galactic nuclei of RE J1034H-396. In 
case of the former object, the source is known because of 
displaying a strong X-ray aperiodic variability with a variety of 
low-frequency (0.08-18 Hz) and high-frequency QPOs (100- 
285 Hz) [TSl-f^S]. Here, we examine its activity peaked at the 
frequency of 4 Hz and derive corresponding QPO significance 
in the wavelet domain. As for the latter source, we select the 
active galaxy of RE J 1034-1-396 which has been recently found 
to uncover the very first ever significant signature (over 3ct) 
of quasi-periodicity among all active galaxies ||49l . Inspired 
by Fourier results of |49|, we verify their findings by determi- 
nation of wavelet significance for detected ~ 2.7 x 10^'* Hz 
variability. 

A. XTE J1550-564 

In order to examine XTE J 1550-564 light curve variability 
which has been dominated by intrinsic red-noise component in 
IG-^-lQi Hz frequency band, after [501 we select RXTE/PCA 
observation of 40099-01-24-01 (1998-09-29) from the public 
HEASARC archive. We reduce the data with LHEASOFT 
package ver 6.6.1 applying the standard PCA selection criteria 
for above data set (see http://heasarc.gsfc.nasa.gov for details). 
As the final product we extract a 2-13 keV light curve with a 
bin size of At = 2^^ s for the first part of PCA observation 
(i.e. N = 80065, T = NAt = 2502 s). 

Fig. [1] presents Fourier power spectrum density calculated 
from 78 averaged spectra based on 1024 point data segments. 
The overall PSD shape can be split into two main compo- 
nents: power-law underlying spectrum and broad line profile 
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Fig. 1. Power spectrum density of XTE J 1550-564 as observed by 
RXTE/PCA on 1998-09-29. Black solid line denotes fitted model composed of 
double broken power-law model (red solid line) and Lorentzian (dotted line) 
representing red-noise background spectrum and quasi-periodic oscillation 
(QPO), respectively. Power spectrum for simulated TK time-series has been 
marked in blue (circle markers). 
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representing red-noise source variability and QPO component, 
respectively. In the fitting process we found that a model 
composed of three-segment broken power-law (i.e. with two 
break frequencies) and Lorentzian profile line 15 1 1 fit the PSD 
best yielding reduced-x^ of 10.9 at 59 degrees of freedom. The 
corresponding power-law slope values and break frequencies 
have been determined as follows; ai — 0.22 ± 0.02, 012 = 
0.90 ± 0.02, ag = 3.40 ± 0.07, /i = 1.63 ± 0.08 Hz and 
/2 = 8.12±0.10, whereas for Lorentzian the central frequency 
of /o = 4.104±0.006 Hz and its FWHM A/ = 0.33±0.01 Hz 
were fitted providing QPO quality factor of (5 — 12.4. We 
marked both model components in Fig. [T] by red solid line 
and dotted line, respectively, where the overall best model has 
been indicated by black solid line. 

Based on the above power-law model we generated a 2^^ 
point long simulated signal (64 segments of 1024 points) with 
TK method. In simulation we assumed the sampling time of 
At — 2^^ s and ensured that resulting rms value (here, derived 
in 0.5-2 Hz and 8-12 Hz frequency band; see Sect. III-BI for 
details) was nearly the same as compared to corresponding rms 
values for X-ray source. We calculated the averaged Fourier 
PSD for TK time-series and plotted it in Fig. [Uby blue circle 
markers. One can clearly notice that TK PSD reproduces well 
the assumed double broken power-law model component. 

Following the recipe given in Section III-CI we aim at 
determination of wavelet confidence at 0.9995 level. From the 
entire XTE J 1550-564 light curve, we select a 3072 point long 
test signal (96 s) for which we compute wavelet significance 
matrix (Eq. |20| |. Upper plot of Fig. [2] presents a dependence 
of 0.9995-th wavelet quantile (Eq. [TsT l in a function of Fourier 
frequency in 1-10 Hz band, both for X-ray data segment 
(black solid line) and TK simulation (red dashed line). It is 
straightforward to note that at the significance of 5 x lO^"' 
the quasi-periodic variability around / ~ 4 Hz is high above 
the mean (simulated) red-noise level. Finally, the bottom map 
of Fig. 12] presents the wavelet power spectrum (Eq. [T]! for 
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Fig. 2. (upper) The dependence of 0.9995-th quantiles in a function 
of frequency calculated for selected 96 s data signal of XTE J 1550-564 
(solid black line) and obtained from wavelet analysis of simulated TK series 
(red dashed line). See Section IIII-AI for details, (bottom) Con'esponding 
wavelet power spectrum displayed with derived significance contours (0.9995 
confidence level). For better clarity of map reading we plot 96 s long wavelet 
spectrum between 37-94 s only. 

considered data segment with overlaid contours of derived 
significance. 

B. REJ1034+396 

A Seyfert 1 galaxy of RE J1034H-396 has been found 
and confirmed as the first active galaxy showing highly 
significant quasi-periodic oscillation (QPO) at the frequency 
/ ~ 2.7 X lO-"* Hz El. We reproduce their PSD in Fig. [3]as 
computed for signal segment of XMM-Newton EPIC observa- 
tion of RE J1034H-396 conducted on 2007-05-31. SoHd black 
line crossing PSD denotes best-fitted power-law spectrum 
described as: 

lgP(/) = -1.35 Ig/- 4.128 (21) 

(M. Gierlinski, private communication). Upper two dashed 
lines mark computed confidence levels for detected period- 
icity at 99.73% (3(7) and 99.99% confidence levels. In order 
to verify aforementioned significance of QPO detection, we 
downloaded from XMM-Newton Science Operations Centre 
the same EPIC pnH-MOSl/MOS2 observation data set and 
reduced it in the same manner as described in Il49l . We limit 
original X-ray light curve (0.3-10 keV) to 601 point long time- 
series (corresponding to segment 2 of 1491 ) where the sampling 
time of A< = 100 s has been used. 
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Fig. 4. RE J 1034+396: Wavelet power spectrum calculated for the second 
data segment of XMM-Newton EPIC/pn+MOS observation, as defined by 1 49 1 . 
Quasi-periodic modulation is clearly detected at / ~ 2.7x 10^* Hz frequency 
and last throughout a whole light curve duration. Outer and inner contours 
correspond to derived 99.73% (3cr) and 99.99% confid ence l evel, respectively, 
with the use of new algorithm as provided in Section Hl-C I 



Fig. 3. Fourier power spectral density of an active galaxy of RE J 1034+396 
|49 1 displaying a significant quasi-periodicity at frequency / ~ 2.7x 10^* Hz 
(> 99.99% of confidence; upper dotted line). Solid line denotes best power- 
law model fitted to the data P{f) oc f~^'^^, whereas bottom horizontal 
dotted fines marks expected white-noise level. Reproduced with kind permis- 
sion of Marek Gierlinski. 

We compute wavelet power spectrum ([1]) in 10~^-10~'^ Hz 
band and make use of our new algorithm in order to determine 
two contour plots marking 99.73% and 99.99% confidence 
for wavelet results. For the purpose of wavelet background 
simulation, we generate 2^^ point TK signal based on assumed 
single power-law model (Eq. 1211 1. Since the underlying PSD 
model is known in advance, we perform an integration of it 
between 10^^ Hz and 10"'^ Hz to find the corresponding 
rms value in this frequency range (rms~ 5.5%). Within the 
simulation process of TK time-series we check to make sure 
that in the same frequency band the rms of TK signal stays 
in agreement with derived value. Based on that, we determine 
two significance matrices (Eq. |20b and plot them together with 
wavelet power spectrum. The final form of all calculations 
presents Fig. |4] where an outter and inner contour denote 
99.73% and 99.99% confidence level, respectively. 

Remarkably, a detected quasi-periodic oscillation at / ~ 
2.7 X 10^^ Hz appears to be very coherent in time though 
showing small but noticeable period drift. We checked that the 
entire periodicity remains significant at the confidence level of 
99% throughout the signal segment duration. In contrast, as 
marked in Fig.|4] only a small fraction of QPO activity exceeds 
top level of 0.01% significance whereas about 70% of detected 
QPO signal stays significant at 3(j level. 

IV. Conclusions 

In this paper we reviewed a recent status of knowledge 
devoted to wavelet significance tests and proposed a new 
method designed especially for signals which intrinsic vari- 
ability displays in Fourier power spectral density a power-law 
form. Since a large number of astronomical X-ray sources has 
been found to be characterized by the similar PSD shapes, we 



focused our attention to that class of objects and their time- 
series analysis. 

We presented a flexible algorithm helpful in determination 
of the significance of quasi-periodic oscillations in the wavelet 
domain. Our method has been based on a comparison of 
wavelet coefficients calculated for X-ray time-series with those 
computed for simulated signal, where the latter displays the 
same PSD broad-band model of variability as X-ray object. 

We argue that the usefulness of our new algorithm can be 
easily extended to non-astronomical time-series which PSD 
can be described by power-law model as well. If required, 
the omission of exponential transformation of linear signal 
generated with the use of TK method and inclusion of Poisson 
noise should be avoided (Section [II- Al l. 
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